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Abstract 

It is known that relative entropy of entanglement for an entangled state p is defined via its 
closest separable (or positive partial transpose) state a. Recently, it has been shown how to find 
p provided that a is given in two-qubit system. In this paper we study on the reverse process-i.e., 
how to find a provided that p is given. It is shown that if p is one of Bell-diagonal, generalized 
Vedral-Plenio, and generalized Horodecki states, one can find a from a geometrical point of view. 
This is possible due to the following two facts: (i) The Bloch vectors of p and a are identical 
with each other (ii) The correlation vector of a can be computed from a crossing point between a 
minimal geometrical object, in which all separable states reside in the presence of Bloch vectors, 
and a straight line, which connects the point corresponding to the correlation vector of p and the 
nearest vertex of the maximal tetrahedron, where all two-qubit states reside. It is shown, however, 
that these nice properties are not maintained for the arbitrary two-qubit states. 
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I. INTRODUCTION 



It is well known that entanglement of quantum states is an important physical resource 
in the context of the quantum information theories. It plays a crucial role in quantum 
teleportation 1|, superdense codingj^J, quantum cloning 3], quantum cryptography 4] and 
quantum computer technology j^, Q] 1 . Therefore, to understand how to quantify and how 
to characterize the entanglement for a given quantum state is an highly important physical 
task. 

Many entanglement measures have been developed for last few years. Above all, in 
our opinion, the most important entanglement measure is a distillable entanglement [lOj, 
which quantifies how many maximally entangled states can be constructed from the copies 
of the given quantum state in the asymptotic region. The importance of the distillable 
entanglement arises due to the fact that the entanglement is fragile when noises interfere 
the quantum information processing. The disadvantage of the distillable entanglement is 
its calculational difficulty In order to compute the distillable entanglement analytically 
we should find the optimal purification (or distillation) protocol. If this optimal protocol 
generates n maximally entangled states from m-copies of the quantum state p, the distillable 
entanglement for p is given by 2 

n 

Dip) = lim — . (1.1) 

m— ¥oo Tfl 

However, finding an optimal purification protocol is an highly non-trivial task. It makes it 
difficult to compute the distillable entanglement analytically. 



Fortunate 
in Ref. 



11 



y, the tight upper bound of the distillable entanglement has been developed 



121 ] . In these references the new entanglement measure called relative entropy of 



entanglement (REE) was introduced. It is defined as 



E R (p)=minS(p\\a), (1.2) 



1 There are, however, several examples, where entanglement does not play an important role in the quantum 
computation. For example, the efficiency of Grover's search algorithm gets worsened if the initial state 
is entangled oneQ. Another important example is a deterministic quantum computation with one pure 
qubitQ. Other simple examples are presented in Ref.Q. Therefore, one cannot conclude definitely that 
entanglement is essential for quantum computation. 

2 In Rcf.[10( the distillable entanglement D is divided into D\ and D2 depending on one-way and two-way 
classical communications. Throughout this paper we only consider the two-way classical communication. 
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where V is a set of separable states and S{p\\a) is a quantum relative entropy-i.e., S(p\\a) 



tr(plnp — plna). It was shown in Re: 



entanglement. Subsequently, Rains 
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121 ] that Er(p) is a upper bound of the distillable 
141 ] has shown that 



E R (p) = min S{p\\a) 



;i.s) 



where T>ppt is a set of positive partial transposition (PPT) states, is more tight upper 
bound when p is an higher-dimensional bipartite state. Using the facts that the REE is 
an upper bound of the distillable entanglement and Smolin state lj]] is a bound entangled 
state, the distillable entanglement for the various Bell-state mixtures has been analytically 
computed KMlq]. In order to understand the distillable entanglement more deeply, therefore, 
it is important to develop the various techniques for the explicit computation of the REE. 
Of course, regardless of the distillable entanglement, the development of the calculation 
technique for the REE itself is important to understand the characterization of entanglement 
more profoundly. For last few years many properties of the REE were investigated [ljj. 
Furthermore, recently relation between the REE and other distance measures has been 



studied 



20 



2l|. 



In this paper we confine ourselves to the REE when p is two-qubit states-i.e., p G 'H 2 ^)'H 2 . 
Since there is no bound-entangled state in this case, E R {p) and E R (p) defined in Eq. fll.2p 
and Eq. (11.31) are same. Let a* be the closest separable state (CSS) of p. Then, E R (p) is 
given by 

E R (p)=wmS(p\\a) = S(p\\o*). (1.4) 
When the CSS cr* is explicitly given and it is full rank, Ref. 22J has presented how to 



construct the set of the entangled states, whose CSS are a* . Let \i) and Aj be eigenvectors 
and corresponding eigenvalues of cr*. If a* is the CSS (hence, edge) state, then its partial 
transposition cr r is rank deficient. Let \<f>) be the kernel of 



a 



0. 



;i.s) 



Then, the set of the entangled states p(x), whose CSS are a*, is given by the following 
one-parameter family expression: 



p{x) = cr* - xG(a*) 

G(a*) = ^G^|.)(.|(|0)(0|) r |j)(j|, 



;i.6) 



where x > and 



G 



'j 



for i = j 
for i 7^ j. 



;i.t) 



In Xi— In Aj 

When, however, the entangled state p is explicitly given, it is difficult to use Eq.f ll.6p for 
finding its CSS. In other words we have to find the reverse process of Ref.jzj] in order to 
derive the closed formula of the REE for the arbitrary two-qubit state p as Wootters .2 3] h as 



done in the entanglement of formation. Unfortunately, still it is an unsolved problem 



In this paper we will explore the reverse process of Ref. 22|. We will show that the 



24] 



reverse process of Ref. 



221 ] is possible, at least, for the Bell-diagonal, generalized Vedral- 



Plenio, and generalized Horodecki states. We present a method for finding the corresponding 



CSS systematically 



Horodecki in Ref. 



or these states by generalizing the geometrical method discussed by 



25]. We also discuss why it is difficult to find the CSS for the arbitrary 



25] such as tetrahedron T and 



two-qubit mixed states from the geometrical point of view. The paper is organized as follows. 
In section II we will show how to find the CSS for the Bell-diagonal states. In section III 
we discuss how the geometrical objects presented in Ref. 
octahedron C are deformed in the presence of the non-zero Bloch vectors. In section IV 
and section V we present a method for finding the corresponding CSS for the generalized 
Vedral-Plenio and generalized Horodecki states, respectively. In section VI we discuss why it 
is very difficult task to find the CSS for the arbitrary two-qubit states from the geometrical 
point of view. In section VII a brief conclusion is given. 



II. CSS FOR THE BELL-DIAGONAL STATES 



In this section we show how to find the CSS when p is the Bell-diagonal state from the 
geometrical point of view. In fact, this problem was already solved in Ref. 12| long ago. The 
reason why we re-consider the same problem is to stress the geometrical analysis. 

An arbitrary two-qubit state can be represented as follows: 

3 

I®I + r- (T(3l + I(3S-(T+ ^ 9r, 



P 



cr r , 



m,n=l 



(2.1) 



where r and s are Bloch vectors and Oi is usual Pauli matrix. The coefficients g mn form a 
real matrix and represent the interaction of the qubits. If state p is explicitly given, one can 
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derive the Bloch vectors r and s and the correlation tensor g^ as follows: 
r = ti(p A cr), s = tr(pBtr), g {j = tx{pui <g> cr,-), 



(2.2) 



where pA = tr^p and ps = tr^p. It is well known that an appropriate local unitary 
(LU) transformation of p can make g mn to be diagonal (see appendix of Ref.j^f]]). Since 
entanglement is invariant under the LU transformation, it is in general sufficient to consider 
the case of diagonal g mn for the discussion of entanglement. Thus, without loss of generality, 
one can express p as 



P = 

If p = where 



I<g>I + r- <T®I + I<g)s-(T + g n a n <S> a r , 



n=l 



V2 
1 



(|00> + |11» 



V2 



(|00>-|11» 



(2.3) 



(2.4) 



| j 9 8 > = -^(|01> + |10» 



|/3 4 ) = _(|01>-|10», 



it is easy to show that the corresponding Bloch vectors r and s are vanishing and the 
corresponding correlation tensor g mn become 

0i = diag(l, -1, 1) g 2 = diag(-l, 1, 1) g 3 = diag(l, 1, -1) g 4 = diag(-l, -1, -1). 

(2.5) 

If, therefore, p is Bell-diagonal state, the Bloch vectors r and s are always null vectors. 

Since we are considering on the diagonal case of the correlation tensor, we will regard, 
from now on, the tensor as a vector, whose components are equal to the diagonal elements. 
When r = s = 0, Horodecki has shown in Ref. 25j that the total two-qubit states belong 
to the tetrahedron T with vertices v\ = (1, —1, 1), y>i = (—1, 1, 1), v 3 = (1, 1, —1) and v& = 
(— 1, —1, —1) in the correlation vector space. Ref. 25] also has shown that the separable states 
(with r = s = 0) belong to the octahedron C with vertices = (±1, 0, 0), = (0, ±1, 0) 



(o[-\o ( f\o ( ^) and 



and 03 = (0, 0, ±1). This is pictorially represented in Fig. 1. 

As Fig. 1 shows, the planes (oi + \ o 2 ■*), {o^\o^\o^) 
(o[ ,02 ,03 "*) are parts of the planes (fi, t>3, 1*4), (i>i, U2, U3), (vi,V2,Va) and (^2,^3,^4) Re- 
spectively 3 . Therefore, all entangled Bell-diagonal mixtures belong to the small four tetra- 
hedra (i*,«4 + \ojf\<4 +) ), («a, <t\ 4 +) , 4 +) ), (v 3 ,o{ +) ,o { 2 +) ,0^) and (v^A^A^)- 



This statement can be confirmed by deriving the respective plane equations. The plane equations for 
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FIG. 1: The total Bell-diagonal states belong to the tetrahedron (1)1,1)2,1)3,1)4) and the set of the 
separable states belong to the octahedron, whose vertices are of, of and of . As this figure shows, 
the planes (oj , o 2 ,o 3 ^ ) , (o^\ o^\ o^ ) , (o[ ,o 2 , 03 ) and (o[ ?°3 ') are contained in 

the planes (vi, v%, V4), (vi,V2,vs), {v\,V2,v^) and (^2,^3,^4), respectively. Therefore, all entangled 
Bell-diagonal mixtures belong to the small four tetrahedra (v\ , o^ ,o 2 \ o 3 + ' ) ) , [1)2 , o\ \ o 2 + ^ , o 3 + ^ ) , 

(^3,0 1 +) ,4 +) '4~ ) ) and ( v 4,o[~\o { 2 ~ ) ,0 ( f ) ). 



Now, we show how to perform the reverse process of Ref. 22j when p is an entangled 
Bell-diagonal state. This can be achieved by following two theorems. 

Theorem 1. Every Bell state has infinite CSS, which cover fully the nearest surface of 
the octahedron L. 

Proof. It is sufficient to prove this theorem when p = When p = \/3i)(/3i\ (i = 

2, 3, 4), one can prove the theorem similarly. 
Let a be a following Bell-diagonal state: 

3 



1 



n=l 



(2.6) 



with p = (x, y, z). Then, it is easy to show that the spectral decomposition of cr is 



a 



4 

1 + x + y 



(2.7) 



X 



y 



WW- 



4 • 



The nearest surface of C from p = I is (or\o : 



(+) „(-) „(+) 



>°3 



whose surface equation is 



x—y+z = 1. If a belongs to the surface (o^\ o 2 , o 3 + ' ) ), it is easy to show that S(p\ \a) = In 2, 



which exactly coincides with the REE of the Bell states Therefore, a on the surface 
( 0l +3 \o\~ ) ,of ) ) is the CSS of |/3i)(/?i|. 

Now, let us consider the case where a belongs to other surface. For example, let us assume 
that a belongs to the surface (o^\ o^, o 3 ^), whose surface equation is x + y — z — 1. Then, 



(vi, V3, V4), (vi,V2, V3), (v 1 ,v 2 ,V4) and (v 2 ,v 3 , V4) are x — y — z = 1, x + y + z = 1, -x — y + z = I and 
—x + y — z = 1, respectively. It is easy to show that these plane equations are the same planes with the 
planes {o{ +) , of\ o^ } ), (o[ +) ,o { 2 +) ,o { 3 +) ), (o[ -) , of\ o| +) ) and (o { ~ ] , o| +) , o^), respectively. 
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S(p\\a) reduces to In 2— In x, which is less than In 2 if x ^ 1. Therefore a on (o[ , c4 , Og"" 1 ) is 
not CSS of By same way one can show that a on (oj , o^, Og + ' ) ) or (oj~ , Oj , ojj ) 

is not CSS of which completes the proof. 

Theorem 2. TTie GS'S' o/ the any Bell- diagonal state p corresponds to the crossing point 
between the nearest surface of C from p and the straight line I, which connects p and the 
nearest vertex ofT from p. 

n 

Proof. If a is CSS of p, the CSS of p = xp+ (1 — x)o~ is also 0"ll2|. Let p be p = 
Then, theorem 1 implies that a can be any point on the surface (o^\o^~\o^). Let p 
belong to the small tetrahedron (v\, o± , \ <4 )■ Note that p corresponds to a internally 
dividing point of the line segment pa. Since Eq. fll.6j) implies that the set of the entangled 
states which have same CSS should be represented by the straight line, the only possible 
a as CSS of p is a crossing point between a line pp and the surface (o[ + \ o^\ o!^), which 
completes the proof for the Bell-diagonal states. 

By making use of the Theorem 2 one can always find the CSS a if p is a Bell-diagonal 
state. Fig. 2 shows how to find the CSS for the Bell-diagonal state. First, extend the 
line segment between p and the point corresponding to the nearest vertex of T. Second, 
compute the coordinate of the crossing point between the line and the nearest surface of the 
octahedron C. Finally, find the CSS which corresponds to the crossing point. This complete 
the reverse process of Ref. 22] . 



III. GEOMETRICAL DEFORMATION OF T AND C 

When the Bloch vectors r and s are non-zero, the geometrical objects T and C should be 
deformed. In this section we will discuss how T and C are deformed. In order to perform the 
following analysis analytically we consider in this paper the case where r and s are parallel 
to each other. It is worthwhile noting that if r and s are x- or y-direction, one can make 
them to be ^-directional via the appropriate local-unitary transformation. For example, if 
they are x-direction, p' = (U <8> U)p{U <g> U)^ with 
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V2 



FIG. 2: This figure shows how to find the CSS for the Bell-diagonal state. First, extend the line 
segment between p and the point corresponding to the nearest vertex of T ■ Second, compute the 
coordinate of the crossing point between the line and the nearest surface of the octahedron C. 
Finally, find the CSS of p which corresponds to the crossing point. 

has ^-directional Bloch vectors and its correlation vector changes from g 2 , g%) to 
{~93i 92, 9i)- Similarly, one can change the state with ^/-directional Bloch vectors into the 
state with ^-directional Bloch vectors without altering the diagonal property of the correla- 
tion term. 

In this reason it is reasonable to assume that the directions of the Bloch vectors are 
z-direction by writing r = (0,0, r) and s = (0,0, s) 4 . In this case the arbitrary two-qubit 

4 Even if r and s are not parallel with each other, one can make them to be z-directional via an appropriate 
local-unitary transformation. In this case, however, the correlation term loses its diagonal property. 
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state p defined in Eq. fl2.3p with g = (q 1} q 2 , q 3 ) reduces to 



1 + r + s + q 3 



?1 - <?2 




1 + r - S - q 3 qi+q 2 
qt + q 2 l-r + s-q 3 



q% ~ <?2 





(3.1) 











1 - r - s + q 3 J 



FIG. 3: The deformation of T is plotted when r = s = 0.3 (Fig. 3a), r = — s = 0.3 (Fig. 3b), 
r = s = 0.5 (Fig. 3c) and r = — s = 0.5 (Fig. 3d). For comparison we plot T together. The 
appearance of non-zero Bloch vectors generally shrink the tetrahedron. The shrinking rate becomes 
larger with increasing the norm of the Bloch vectors. 



The eigenvalues and eigenvectors of p are summarized at Table I. 



eigenvalues of p 


eigenvectors of p 


p± = \{(l-q 3 )±M l } 
v± = \{{l + q 3 )±M 2 } 


K> = 


i 

A± 
1 
A± 


[(gi + ? 2 )|01>-{(r-s) 
[(9i-(fe)|00)-{(r + s) 


TMi}|10)] 
TM 2 }|11)] 



Table I: Eigenvalues and Eigenvectors of p in Eq. fl3.ll) 



At Table I Mi, M 2 , A± and A± are given by 

Mi = ^( r - s ) 2 + (qi + q 2 y M 2 = y/(r + s) 2 + ( Ql - q 2 f (3.2) 
A ± = ^/{(r-s)TM 1 } 2 + (q 1 + q 2 y A± = ^{(r + 7) T M 2 } 2 + ( 9l - q 2 f. 

Then, the deformation of T can be obtained from the positivity condition of p. Since 
deformation should be a set of the boundary states, the condition of the deformation becomes 

min(//_, V-) = 0. (3.3) 

One can make two surfaces by making use of Eq. (13.31) . Each surface corresponds to 
min( / u_,z/_) = /i_ = or min(/i_,z/_) = V- = 0. Gluing these surfaces together yields 
the deformation of T. 

In Fig. 3 we plot the deformation of T when r = s = 0.3 (Fig. 3a), r = —s = 0.3 (Fig. 
3b), r = s = 0.5 (Fig. 3c) and r = — s = 0.5 (Fig. 3d). For comparison we plot T together. 
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For convenience, we will call the deformation of T with fixed r and s as %, s - From Fig. 3 
one can realize that the deformation % )S has following two characteristics. First one is that 
the effect of the non-zero Bloch vectors is to shrink the geometrical object. The shrinking 
rate becomes larger with increasing |r| and \s\. When r = s, the deformation is biased 
toward (t>i,t>2) region. When, however, r = — s, the deformation is biased toward (f 3 ,f 4 ) 
region. The shrinkage of % yS implies that the number of proper quantum states reduces 
with increasing |r| and \s\ due to the constraint trp 2 < 1. The second characteristic of % )S 
is that it has continuous smooth surface while T has sharp edges. This fact arises from 
the condition min(/i_, vJ) = 0. When r = s = 0, this condition generates the four surface 
equations 

± (<?i + 92 ) + <7s = 1 ± (<?i - 92) - q 3 = 1, (3.4) 

each of which corresponds to the surface of T. When, however, r and s are non-zero, these 
four equations reduce to the following two equations: 

y/(r - s) 2 + ( qi + q 2 Y + q 3 = 1 a/ (r + s) 2 + (qi — q 2 ) 2 — g 3 = 1. (3.5) 

This implies that the deformation % tS can be formed by attaching two smooth surfaces when 
r ^ ±s. 

FIG. 4: The deformation of C is plotted when r = s = 0.3 (Fig. 4a), r = — s = 0.3 (Fig. 
4b), r = s = 0.5 (Fig. 4c) and r = —s = 0.5 (Fig. 4d). For comparison we plot C together. 
The appearance of non-zero Bloch vectors generally shrinks the octahedron. The shrinking rate 
becomes larger with increasing the norm of the Bloch vectors 

Now, we discuss the deformation of C when the Bloch vectors are r = (0,0, r) and 
s = (0,0, s). We will call this deformation as C TyS . We assume that p in Eq.f l3.ip is a 
separable state. In this case the PPT state of p, say p r , should be positive. The eigenvalues 
and the corresponding eigenvectors of p r are summarized in Table II. 



eigenvalues of p T 


eigenvectors of p T 


p v ± 


= \{(l-q*)±M*} 


1/4) = 


^[( gi - ft )|01>-{(r- a ) 


T-MfjllO)] 


-I 


= H(l + <? 3 )±M 2 r } 


Wl) = 


^ [(gi + fe)|00)-{(r + s) 


T-M 2 r }|ll>] 



Table II: Eigenvalues and Eigenvectors of p r 
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At Table II M[, M 2 r , A£ and A£ are defined as 



Ml = Mi 



Ml = M 2 



<?2^-<?2 



V(r - s) 2 + (ft - ft) 2 (3.6) 
V(r + s) 2 + (ft + ft) 2 



A r = ^ {(r — s) zp Mf} 2 + ( ft - ft) 2 A£ = ^{(r + s) =F M 2 r } 2 + (ft + ft) 2 . 

Then £ r s can be obtained from the positivity condition of p r . Since, furthermore, £ rjS 
should be a set of the edge states, the condition for the deformation of £ reduces to 

min(//,z/) = 0. (3.7) 

As deformation of T Eq. fl3.7j) generates two surfaces, each of which corresponds to 
min(/z£,z/£) = pF_ = or min = z/ r = 0. Gluing these two surfaces one can 

make the deformation of C. 

In Fig. 4 we plot the deformation of £ at r = s = 0.3 (Fig. 4a), r = — s = 0.3 (Fig. 4b), 
r = s = 0.5 (Fig. 4c) and r = —s = 0.5 (Fig. 4d). For comparison we plot £ together. 
Like the deformation of T, Fig. 4 indicates that the effect of the non-zero Bloch vectors 
is to shrink £ toward a particular direction. Fig. 4 also shows that the shrinking rate 
becomes larger and larger with increasing the norm of the Bloch vectors. Like again, 
the deformation of the octahedron £ rs also has smooth surfaces while £ has sharp edges. 



IV. CSS FOR THE GENERALIZED VEDRAL-PLENIO STATES 

In this section we show how to derive the CSS for the Vedral-Plenio (VP) states. The VP 
states are defined as mixture of one Bell state and separable states, which are not orthogonal 
to the Bell state. One but most general example of the VP state is 

Pv P = A 1 |/3 1 )(/3 1 |+A 2 |00)(00|+A 3 |ll)(ll|, (4.1) 

where = (l/>/2)(|00> + |11» and Ai + A 2 + A 3 = 1. 
Let the arbitrary VP state be 

Pv P = ^ M ® / + r vp ■ cr <g> I + I <g> s vp ■ cr + y^Xt vp ) n a n <g) cr„J (4.2) 



11 



and its CSS be 

n vp = i ^7 <g> 7 + u vp ■ a <g> 7 + 7 <g> U up ■ cr + y^(r„ p ) n cr w <g) cr^j . (4.3) 

The following theorem shows how to compute u vp , v vp and T vp from p vp . 

Theorem 3. If 7t vp is the CSS of p vp , u vp = r vp and v vp = s vp . Let £ be a straight line, 
which connects t vp and the nearest vertex of T . Then, r vp is a crossing point between £ and 
r 

Proof. We will prove this theorem as following procedure. First, we assume that this 
theorem is correct. Then, following this theorem one can derive the trial CSS state of p vp . 
Next, by making use of Eq. (11.61) we will show that this trial CSS state is a really CSS state 
of p vp . 

Since other VP states can be derived from Eq. ( 14. ip by local- unitary (LU) transformation, 
it is sufficient to show that the CSS of p vp in Eq. (14.1 ft satisfies this theorem. The other 
case can be proven similarly. For p vp in Eq.l l4.ip r vp , s vp and t vp become r vp = (0,0, r), 
s v P = (0, 0, s) and t vp = (t 1 , t 2 , t 3 ) where 

r = s = A 2 - A 3 ti = -t 2 = \i t 3 = 1. (4.4) 

Then, a point P = (?i, 92,53) on the line £ satisfies g 2 = —Qi and ^3 = 1. If the point 
P = (91,52,93) is a crossing point between £ and £ riS , the corresponding separable state 
satisfies 

^- = ~l9i|, z/ r = i(l-|A 2 -A 3 |), (4.5) 

where pF_ and v v _ are defined at Table II. Therefore, the CSS condition ( 13 .TP implies q\ = 0, 
which results in T vp = (0, 0, 1). If, therefore, this theorem is correct, the CSS of p vp is 

Uvp = \ [7 <g> 7 + (A 2 - A 3 ) (a 3 ® I + 7 ® a 3 ) + a 3 ® a 3 ] (4.6) 

/^ + A 2 00 \ 



V oof + a 3 / 

In order to show that ir vp in Eq. (l4.6p is really CSS of p vp , it is convenient to define another 
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edge state 



7T 



vp 



(t ® a x )n vp (t <g> a x 



it 



f + A 2 







e 




f + A 3 



(4.7) 





\0 e / 

where the infinitesimal positive parameter e is introduced for convenience. This parameter 
will be taken to be zero after calculation. 
Let us define a edge state 

/ R x \ 

R 2 Y 

y R 3 

\ RiJ 



<5Z 



(4.8) 



with y = y/R~iR~l and -R2-R3 > R1R4. Then, by making use of Eq. fll.6p Ref. 22] has shown 
that the set of the entangled states, which have <Jz as CSS, is represented as 

( Rx-xRi \ 

R 2 -xR 2 Y-xY 

Y-xY R 3 -xR 3 

Ri-xRiJ 



Pz(x) 



\ 








(4.9) 



where x > and 5 



y 2 



R\ + i?4 

-R3 = — 2i?i — i?2 



i? 2 = 2y 2 rf [(R 2 - R 3 )(R 2 L -z) + 2Y 2 L] 
Y = Yd [2Y 2 (R 2 + R 3 )L + (R 2 - R 3 ) 2 z] . 



(4.10) 



In Eq.f ETTO]) we define 

z= ^(R 2 -R 3 ) 2 + 4R!R 4 L 



In 



(R 2 + R 3 ) + z 
(R 2 + R 3 )-z 



d 



1 



(4.11) 



Now, we identify a z with n vp by putting Ri = i? 4 = Y 



(Rx + R^z 2 ^ 
e, R 2 = Ai/2 + A 2 and 



R3 = Ai/2 + A3. Then, it is straightforward to compute R\, R 2 , R 3 , R4 and Y. After taking 
e — > limit, one can show Ri = R 2 = R 3 = R 4 = and Y = — \\ 2 — \ 3 \/2L, where 



(4.12) 



We corrected the sign mistake of Ref. [25 
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Therefore, the set of the entangled states, which have ir vp as CSS, can be represented by 



/o 



Pz{x) 



| + A 2 
x^ ^ + A 3 



\0 



(4.13) 



0/ 



Finally, the set of the entangled states, which have n vp as CSS, can be derived by taking the 
inverse LU transformation, i.e. 



Pz {x) = (1 ® (7 a; ) t / Z (x)(l g> (Xj.) . 



(4.14) 



It is easy to show that pz{x) reduces to p vp in Eq.f l4.ip when x = x vp = A1-L/IA2 — A3 1 > 0, 
which completes the proof. 



FIG. 5: Fig. 5 shows how to find the CSS for the VP states. Fig. 5a and Fig. 5b correspond to 
r = s = 0.3 and r = s = 0.5, respectively. To find a CSS make a straight line I first, which connects 
the nearest vertex of T and a point t vp . Second, compute the coordinate for the intersection point 
between I and C TjS . Thirdly, identify the crossing point with r vp . Keeping u vp = r vp and v vp = s vp , 



one can find the CSS of the VP state. 



Fig. 5 shows how to find the CSS for the VP state geometrically when r = s = 0.3 (Fig. 
5a) and r = s = 0.5 (Fig. 5b). Fig. 5 indicates that the generalized VP states are on the 
edges of T . First we make a line, which connects the nearest vertex of T and p vp . Then, we 
compute the coordinate of the crossing point r between the line and £ r>s . Finally, the CSS 
ii vp of p vp can be computed by Eq.f l4.3p with keeping the Bloch vectors. 

V. CSS FOR THE GENERALIZED HORODECKI STATES 

In this section we discuss how to derive the CSS of the generalized Horodecki states. The 
Horodecki states are defined as mixture of one Bell state and separable states, which are 
orthogonal to the Bell state. One but most general example of the VP state is 



Ph = A 1 |/3 1 )(/3i|+A 2 |01)(01|+A 3 |10}(10| 
14 



(5.1) 



with A1 + A2 + A3 = 1. By contrast with the VP state Horodecki state ( 15. ip is separable when 
X\ < 4A2A3. This can be easily understood by computing the concurrence of pn, which is 




, , - v , if Ai > 2^X^ 

C(ph) = < " " 1 " \1± \ (5.2) 

if Ai < 2^\ z . 

Thus, C(ph) becomes zero when A^ < 4A2A3, which indicates that pu is separable in this 
region. 

Let the arbitrary Horodecki state be 



1 I 3 
p H = j \ l®I + r H -cr®I + I®s H -(T + ^2(t H 

\ n=l 



n a n g) cr„ (5.3) 



and its CSS be 



(5.4) 



The following theorem shows how to compute Uh, Vjj and th from p#. 

Theorem 4. I/tth is a CSS of pn, Uh = th andvn = sh- Let I be a straight line, which 
connects t H and the nearest vertex of T . Then, Th is the nearest crossing point between t 
and C TyS . 

Proof. We will prove this theorem by following the same procedure of Theorem 3. 
Since other Horodecki states can be derived from pn in Eq. (15. 11) by LU transformation, it 
is sufficient to show that the CSS of Eq. (l5.1l) satisfies this theorem. By identifying Eq. (I5.1j) 
with Eq. ()5.3p one can easily show that r H , s H and i# become r H = (0,0, r), s H = (0,0, s) 
and tff = (ti,t2,ts), where 

r = -s = A 2 - A 3 h = -t 2 = Ai t 3 = 2Ai - 1. (5.5) 

Then, a point P(qi, qi, 93) on the line £ satisfies q<i = q% and q% = 2q\ — 1. 

Let the point P be crossing point between t and £ r>s . Then, pF_ and v_ for the state 
corresponding to the point P are given by 

1 



2 



\/9i 2 + (A 2 - A 3 )< 



rS = |. (5.6) 



Therefore, the CSS condition min(/£, ify = gives two solutions Pi = (51,-91,251 — 1) 
and P2 = (0, 0, —1), where 

9i = ^(Ai + 2A 2 )(A 1 + 2A 3 ). (5.7) 
15 



Since we have to choose the nearest point from tjj, the solution we want is the former. 
Therefore, th becomes (gi, — gi, 2gi — 1). Then theorem 4 claims that the CSS of pu is 



I ® I + (\ 2 - \ 3 ) {a 3 ® I - I ®a 3 } (5.1 

+giO"i ® e\ — <?io 2 ® o"2 + (2g - 1)0-3 ® o 3 

/(A 1 + 2A 2 )(A 1 + 2A 3 ) (Ai + 2A 2 )(Ai + 2A 3 ) \ 

1 (Ax + 2A 2 ) 2 

4 (Ai + 2A 3 ) 2 

\(A 1 + 2A 2 )(A 1 + 2A 3 ) (\ x + 2A 2 )(Ai + 2A 3 ) / 

Now, we will show that hh in Eq.( l5.8j) is really CSS of pn by making use of Eq. (jl.6p . In 
order to show this we define tth = (1® 0^)^(1 ®a x y. Then by making use of Eq. fl4.8p and 
Eq.f l4.9p it is straightforward to find a set of the entangled quantum states p(x), whose CSS 
are tth- After taking the inverse LU transformation one can derive p(x) = (1 ®a x yp(x) (1 <g> 
a x ). The expression of p(x) is 

/ Y + xr] Y + xr]\ 

Rx-xrj 
R 4 -xt] 



p(x) 



\ Y + XT] 







(5.9) 



Y + xt] J 



where x > and 



R 1 = -(X 1 + 2X 2 ) 2 



1 

4 ( 



R 4 = ^(Ax + 2A 3 ) 2 



(5.10) 



Y 



-(A 1 + 2A 2 )(A 1 + 2A 3 ) r? 



F 2 



When 



;r 



1 /Ai 
V V2 



(5.11) 



p(x) reduces to p# in Eq. fl5.ip . It is easy to prove that xh > if A 2 > 4A 2 A 3 , which is an 
entangled condition for pn- Therefore, theorem 4 is completely proved. 

Fig. 6 shows how to find the CSS for the generalized Horodecki state pn when r = — s = 
0.3 (Fig. 6a) and r = —s = 0.5 (Fig. 6b). If pu is explicitly given, compute Vh, $h and 
tjj- Then make a straight line which connects a point i# and the nearest vertex of T. Find 
the crossing points between the line and L rs . As Fig. 3 shows, there are two intersection 
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FIG. 6: Fig. 6 shows how to find the CSS for the generalized Horodecki states. Fig. 6a and Fig. 
6b correspond to r = —s = 0.3 and r = —s = 0.5 respectively. In order to find CSS one makes a 
straight line I first, which connects the nearest vertex of T and a point tu- Secondly, we compute 
the coordinate for the nearest crossing point between t and £ r ,s- Thirdly, we identify the crossing 
point with th- Keeping uh = fH and vh = sh, one can find the CSS of the Horodecki state. 



points Pi and P2 for the Horodecki states. This is why the CSS condition min(/i^, = 
gives two different solutions. Using the nearest crossing point (Pi in Fig. 6) one can derive 
t h straightforwardly. Finally using Eq. fl5.4p with imposing u H = r H and v H = s H , one can 
derive n H , the CSS of pn- 



VI. DIFFICULTIES IN FINDING CSS FOR ARBITRARY STATES 

In the previous sections we have shown how to find the CSS for the Bell-diagonal, gener- 
alized VP, and generalized Horodecki states. In fact, it is possible to find the CSS because 
those states exhibit the following nice features. Let r, s and t be Bloch and correlation vec- 
tors of those states. Let u, v and r be Bloch and correlation vectors for the corresponding 
CSS of those states. Then, the features are: 

(i) u — r and v = s. 

(ii) t can be computed from the crossing point between the straight line I and the surface 
for a set of the separable states £ r>s . 

However, such simple but nice features are not maintained for the general mixtures. For 
example, let us consider the comparatively simple model introduced in Eq. fl4.8p and Eq. fl4.9p . 
It is straightforward to show that the first property, i.e. u = r and v = s, is not maintained 
unless R2 = -R3 6 . In order to find the CSS for the arbitrary states, therefore, we have to find 
the explicit relations between (r, s) and (u, v). As far as we know, still this is an unsolved 
problem. 

In addition, one can show that the second property is not maintained too for the general 



When i? 2 = R3, one can show [pz(x),az] — 0. Therefore, Ep_(pz{x) <g> p) = Er(pz{x)) + Er{p) for all 



two-qubit mixture p 14|. 
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FIG. 7: Fig. 7 shows that the property (ii) is not maintained for the arbitrary two-qubit states. 
Four figures in Fig. 7 correspond to r = s = 0.3 (Fig. 7a), r = —s = 0.3 (Fig. 7b), r = s = 0.5 
(Fig. 7c) and r = — s = 0.5 (Fig. 7d). For convenience, we plot T and C T)S together in each figure. 
Each line in Fig. 7 represents a set of t, whose CSS has same r. As Fig. 7 has exhibited, not all 
lines pass one of the vertices of T ■ This fact indicates that the nice property (ii) is not maintained 
for the general mixtures. 

mixtures. Using Eq.f ll.6j) we plot the r (correlation vector of CSS)-dependence of t (correla- 
tion vector of the entangled state) with varying the parameter x when min(/z^, 
Since similar behavior arises when min(//£,^£) = [i T _ = 0, we have not included this case 
in Fig. 7. Four figures in Fig. 7 correspond to, respectively, r = s = 0.3 (Fig. 7a), 
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r = —s = 0.3 (Fig. 7b), r = s = 0.5 (Fig. 7c) and r = — s = 0.5 (Fig. 7d). For convenience, 
we plot T and C r ^ s together in each figure. Each line in Fig. 7 represents a set of t, whose 
CSS has same r. As Fig. 7 exhibits, not all lines do pass one of the vertices of T . This fact 
indicates that unfortunately the nice property (ii) is not maintained for the arbitrary states. 

The non-maintenance of the property (ii) can be proved on the analytical ground by 
making use of the simpler model. Let us consider oz in Eq. (I4.8|) and pz{x) in Eq. fl4.9p . Then, 
the Bloch vectors r, s and the correlation vector t of pz(x) are r = (0,0, r), s = (0,0, s) 
and t = (ti,t 2 ,t 3 ), where 

r = (R 1 + R 2 -R 3 - R A ) - x(R 2 - R 3 ) (6.1) 

s = {R 1 -R 2 + R 3 - i? 4 ) + x{R 2 - R 3 ) 

t 1= t 2 = 2Y - 2xY t 3 = (#! -R 2 -R 3 + R 4 ) - AxR x . 

Of course, if we take x — > limit in Eq. (I6.1j) . the corresponding quantities are the Bloch 
vectors and the correlation vector of oz- Now, let us consider another state irz, which 
can be obtained from oz by changing Y — > Y' and R4 — > R[(i — 1, • • • , 4). In order to 



ensure that itz is CSS we require Y' = a/ R[R' a and R' 2 R' 3 > R'iR' 4 . Then, the set of the 
entangled states £,z(x'), whose CSS are n z , can be obtained from pz{x) by changing Y — > Y', 
Ri — > R[(i = 1, • • • ,4) and x — > x' . Thus, Bloch vectors r', s' and correlation vector i! of 
iz{x') are r' = (0,0, r'), s' = (0, 0, s') and t' = {t[,t' 2 ,t' 3 ), where 

r' = (R[ + R' 2 -R' 3 - R' 4 ) - x\R' 2 - R' 3 ) (6.2) 

s' = (R[ -R' 2 + R 3 - R' 4 ) + x'(R 2 - R 3 ) 

t[ =t' 2 = 2Y' - 2x'Y' t' 3 = (R[ -R' 2 -R' 3 + R' 4 ) - Ax'R[. 

Then, it is straightforward to show that the condition t = t' imposes 

Y'(f - f') - 4(y - F')^i / Y(? - f') - 4(Y - Y')Rt 



x = . x 



A(Y'R 1 -YR' 1 ) " A(Y'R 1 -YR' 1 ) ^'^ 

where f = Ri — R 2 — R 3 + R4 and r' = R[ — R' 2 — R' 3 + R' A . Thus, one can compute the 
crossing point t = t' = p 2 , p 3 ), where /ij becomes 

_ 4(yy'^ - Y'YR[) - YY'{f - f ) 

/il " ^ 2 " 2{Y'R 1 - YR[) (6 ' 4) 

4(y - y)^i^i - (?YR'i - f'Y'Rx) 
1X3 ~ Y'Ri - YR[ ' 
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As a special case we consider the Bell-diagonal case by letting Ri = R 4 = Y = 2Ri = 
2R4 = — 2i?2 = — 2i?3 = a, i?2 = R3 = — 2V = b, R'i = R^ = Y' = 2R'i = 27?^ = — 2i? 2 = 
—2R' 3 = a' and R' 2 = R' 3 = —2Y' = b' . Of course, one can show directly that pz{%) and £z{%') 
are really Bell-diagonal states. Using the normalization conditions 2(a+b) = 2(a'+b') = 1, it 
is easy to verify that the crossing point (pi, /i 2 , ^3) is simply /ix = p 2 = 1 and /i 3 = — 1, which 
is one of the vertices of T. It is worthwhile noting that the crossing point is independent of 
particular choice of Bell-diagonal states pz(x) and £z(a/). This fact implies that all straight 
lines, which connect r and t, pass one of the vertices of T, which is consistent with theorem 
2. 

However for the arbitrary mixtures Eq. fl6.4p implies that the crossing point (/ii, p 2 , ^3) 
is dependent on the choice of the entangled states pz{x) and £z(V). This is why the nice 
property (ii), which holds for the Bell-diagonal, generalized VP, and generalized Horodecki 
states, does not hold for the arbitrary mixture as Fig. 7 has indicated. In order to, there- 
fore, derive the closed formula of Er(p) for the arbitrary two-qubit mixtures p, we have to 
understand how the property (ii) is modified when r, s and t are arbitrary. Unfortunately, 
still this is an unsolved problem too. 

VII. CONCLUSION 

In this paper we have considered how to find the CSS in the two-qubit system from the 
geometrical point of view. Of course, one can straightforwardly compute the REE of the 
state p if its CSS is found. Therefore, it is important to develop a technique for finding CSS 
to overcome the calculational difficulty of the REE. Since, furthermore, the REE is a tight 
upper bound of the distillable entanglement, finding CSS is also important to understand 
the nature of the optimal (or near-optimal) purification protocols. 

If p is one of Bell-diagonal, generalized VP, and generalized Horodecki states, we have 
shown how to find the CSS of p, say a, systematically by proving the following nice two 
properties: (i) The Bloch vectors of a are identical with those of p. (ii) The correlation 
vector of a exactly corresponds to the crossing point between the line I and the geometrical 
object C T)S . Using these two properties it is straightforward to find the CSS of p. 

As we have shown in the previous section, however, these two nice properties are not 
maintained for the general two-qubit states. Therefore, in order to derive the closed formula 
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of Er(p) for the arbitrary mixture p we have to understand how these two properties are 
modified when Bloch and correlation vectors of p are arbitrary. The research into these 
issues is in progress and will be reported elsewhere. 

Another interesting issue, which we will go further, is to explore the REE and the dis- 
tillable entanglement for the higher-qubit or qudit systems. Few years ago, the analytical 
expressions of the distillable entanglement are obtained for some higher- dimensional bipar- 



tite states 



16Hl8|. Authors in those references used the upper-bound criterion D < Er and 



separability property of the Smolin's unlockable state[15| in various cuts. We would like to 
modify Eq. fll.6p to be applicable not only for low-rank a* but also for higher- dimensional 
system. If the generalization of Eq. fll.6p is possible, we can use it to compute the REE and 
the distillable entanglement for many more higher-dimensional states. It may enable us to 
understand the nature of the optimal (or near-optimal) purification protocols. This work is 
in progress too. 
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